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Abstract 

We investigate the possibility that Bose-Einstein condensates (BECs) , loaded 
on a 2D optical lattice, undergo - at finite temperature - a Berezinskii- 
Kosterlitz-Thouless (BKT) transition. We show that - in an experimentally 
attainable range of parameters - a planar lattice of BECs is described by 
the XY model at finite temperature. We demonstrate that the interference 
pattern of the expanding condensates provides the experimental signature 
of the BKT transition by showing that, near the critical temperature, the 
= component of the momentum distribution and the central peak of the 
atomic density profile sharply decrease. The finite-temperature transition for 
a "iD optical lattice is also discussed, and the analogies with superconducting 
Josephson junction networks are stressed through the text. 
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I. INTRODUCTION 



Recent advances in the manipulation of cold atomic gases allow nowadays for storing 
Bose-Einstein condensates (BECs) in ID [1-5], 2D [6] and 3D [7] optical lattices, opening 
a pathway for the experimental and theoretical analysis of finite temperature transitions in 
atomic systems. In a way similar to the spectacular realizations obtained with supercon- 
ducting Josephson networks [8,9], arrays of BECs can provide new experimentally realizable 
systems on which to test well known paradigms of the statistical mechanics: our interest here 
is to evidence how to detect a finite-temperature defect-mediated transition in 2D optical 
lattice of BECs. 

According to the well-known Mermin- Wagner theorem [10], two-dimensional systems 
with a continuous symmetry cannot sustain long-range order in the thermodynamic limit at 
finite temperature; however, a phase transition is still possible and it occurs via the unbinding 
of point defects like dislocations or vortices. Defect-mediated transitions are widely studied 
[10] since they provide crucial insights for a wide variety of experiments on thin films. The 
BKT transition is the paradigmatic example of a defect-mediated transition [10] which is 
exhibited by the two-dimensional XY model [11,12], describing 2-components spins on a 
two-dimensional lattice. In the low-temperature phase, characterized by the presence of 
bound vortex-antivortex pairs, the spatial correlations exhibit a power-law decay; above 
a critical temperature Tbkt, the decay is exponential and there is a proliferation of free 
vortices. The BKT transition has been observed in superconducting Josephson arrays [9] 
and its predictions are well verified in the measurements of the superfiuid density in '^He 
films [13]. 

In finite magnetic systems the BKT transition point is signaled by the drop to of a 
suitably defined magnetization [14] ; the existence of such a magnetization in finite systems 
does not contradict the Mermin- Wagner theorem, since the latter is valid only in the thermo- 
dynamic limit. In a large class of experimental situations, the finite size XY model predicts 
magnetization exponents in agreement with the measurements carried out for layered mag- 



2 



nets with planar spin symmetry [15]. We shall show in the sequel that the Fourier transform 

— * 

at /c = of the condensates wavefunction provide the analogous of this magnetization for 
atomic systems, which are naturally finite (since they are trapped by an harmonic potential). 

In this paper we show that - at a critical temperature Tbkt lower than the temperature 
Tbec at which condensation in each well occurs - 2D lattices of BECs may undergo a 
phase transition to a superfluid regime where the phases of the single-well condensates are 
coherently aligned allowing for the observation of a Berezinskii-Kosterlitz-Thouless (BKT) 
transition. We show in the following that the recently realized 2D optical lattice of Bose- 
Einstein condensates [6] may allow for the observation of a BKT transition in a finite bosonic 
system, since the thermodynamic properties of the bosonic lattice at finite temperature may 
be well described, in a suitable and experimentally attainable range of parameters, by the 
Hamiltonian of the XY model. It is easy to convince oneself that the transition we propose 
to observe is different from the quantum phase transition reported in [7] , where the system is 
at T = and the insulator phase (signaled by a reversible destruction of the phase coherence 
across the lattice) is reached varying the optical lattice parameters: at variance, here, we 
propose to fix the system parameters in the superfiuid phase and increase the temperature 
T until the thermally induced vortex proliferation determines the BKT transition. 

We shall also show that the experimental signature of the BKT transition in bosonic 
planar lattices is obtained by measuring, as a function of the temperature, the central peak 
of the interference pattern of the expanding condensates released from the trap when the 
confining potentials are switched off. In fact, the peak of the momentum distribution at 

— * 

A; = may be regarded as the magnetization of a finite size 2D XY magnet, and, as shown 
in the following, it should exhibit a sharp decrease around Tbkt- As a consequence, also the 
central peak of the atomic density profile decreases around Tbkt, signaling the occurrence 
of the BKT transition. 

The plan of the paper is the following: in Section II we discuss the Hamiltonian describing 
a 2D lattice of Bose-Einstein condensates at finite temperature and we show that - in a 
suitable range of parameters - the system is described by the XY Hamiltonian. In Section 
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Ill we show how the signature of the BKT transition may be evidenced in a bosonic planar 
lattice; there we shall also discuss interesting analogies with superconducting Josephson 
junction networks. Section IV is devoted to our concluding remarks. In Appendix A we 
report the details of a variational estimate of the coefficients of the Hamiltonian describing 
ID, 2D and 3D optical lattices, while Appendix B contains the analj^ical computation of 
the Fourier transform of a vortex on a lattice. 

II. THE FINITE-TEMPERATURE HAMILTONIAN FOR THE SYSTEM 

For 2D optical lattices, when the polarization vectors of the two standing wave laser 
fields are orthogonal, the resulting periodic potential for the atoms is 

Voptir) = Vo[sm^ {kx) + sin^ (%)] (1) 

where A; = 27r/A is the wavevector of the laser beams. The potential maximum of a single 
standing wave, Vq = sEr, may be controlled by changing the intensity of the optical lattice 
beams, and it is conveniently measured in units of the recoil energy Er — h'^k'^/2m (m is 
the atomic mass), while, typically, s can vary from up to 30 (gravity is assumed to act 
along the ^-axis). 

Usually, superimposed to the optical potential, there is an harmonic magnetic potential 

Kn(r1 = ^k^(x^ + y^)+c.,V], (2) 

where cu^ {uJr) is the axial (radial) trapping frequency. The minima of the 2D periodic 
potential (1) are located at the points j — (ji,j2) • f with ji and j2 integers, and the 
potential around the minima is Vopt ~ Y^rii^ ~ ^3i/'^Y + {u ~ ^h/'^Y] with 

UJr = ^2Vok^/m (3) 

providing the frequency of the wells. When ujr ujr,u!z, the system realizes a square array 
of tubes, i.e. an array of harmonic traps elongated along the z-axis [6]. Although the 
axial dynamics is very interesting in its own right [16], in the following we analyze only 
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the situation in which the axial degrees of freedom are frozen out, which is reahzable if cuz 
is sufficiently large. We also assume that the harmonic oscillator length of the magnetic 



potential \Jfi/muJr is larger than the size L (in lattice units) of the optical lattice in order to 
reduce density inhomogeneity effects and to allow for a safe control of finite size effects. In 
addition, when all the relevant energy scales are small compared to the excitation energies 
one can safely expand the field operator [17] as 

^(^t)=E ^7^*7(^1 (4) 

3 

with $j(r) the Wannier wavefunction localized at the j'-th well (normalized to 1) and 
N:' — iplip:' the bosonic number operator. Substituting the expansion in the full quantum 

J 3 

Hamiltonian describing the bosonic system, leads to the Bose-Hubbard model [18,17] 

<^,3> 1 

where Z^<j7> denotes a sum over nearest neighbours, U — {Anh^a/m) J df^i (a is the s-wave 



<i,j> 

scattering length) and X ~ - / df[|^ ■ V$j + %Ka:t*jj 



Upon defining 

J = 2KNo (6) 

(where A^o is the average number of atoms per site), when A^o ^ 1 and J/Nq <^ [/, the 
Bose-Hubbard model reduces to 

u 

H^HxY--Y.Q^j (7) 
j 

which describes the so-called quantum phase model (see e.g. [19,20]): in Eq.(7) 9-^ is the 
phase of the j-th condensates and 

HxY^-J E cos(%-^.) (8) 

<h3> 

stands for the Hamiltonian of the classical XY model. When J ^ U and at temperatures 
T 3> U/kB, the pertinent partition function describing the thermodynamic behaviour of the 
BECs stored in an optical lattice may be computed using the classical XY model (8). 



The close relationship of the Bose-Hubbard model (5) and the quantum phase model (7) 
is by now well known [21] and it has been widely used in the analysis of superconducting 
Josephson networks; since, in superconducting networks, Nq is the average number of excess 
Cooper pairs at a given site, the requirement A^o ^ 1 is of course satisfied and the condition 
J/A^Q ^ U may also be easily realized for reasonable values of J. At variance, in bosonic 
lattices A^o varies usually between ~ 1 and ~ 10^ and the validity of the mapping between 
the Bose-Hubbard model and the quantum phase Hamiltonian is not always guaranteed; 
furthermore, in order to get the XY model (8) from the Bose-Hubbard model (5) U should 
be much smaller than J, but not vanishing, in order to satisfy to the condition J/Nq <^ U. 
Of course the Bose-Hubbard model - for [/ = - describes harmonic oscillators and thus 
cannot sustain any BKT transition; nevertheless, when A^o ^ 1 and J/N^ <^ U <S J, 
the Bose-Hubbard Hamiltonian reduces in the XY model (8), which do display the BKT 
transition. 

A simple estimate of the coefficients J and U may be obtained by approximating the 
Wannier functions in Eq. (4) with gaussians, whose widths are determined variationally. The 
details of the computation are reported in Appendix A. In Fig.l we plot J/U as a function 
of Vq for a 2D optical lattice; for comparison, we plot the same ratio for ID and 3D optical 
lattices also. Since, in a mean-field approach, the Mott insulator-superfluid transition occurs 
at T = approximately at 2zJ/U ~ 1 [18,20] {z is the number of nearest-neighbours), 
from Fig.l one sees that it is much easier to detect a quantum phase transition for the 
three-dimensional array (and indeed it has been recently detected [7]); no quantum phase 
transition has yet been reported in literature for 2D and ID optical lattices, since a larger 
laser power Vq is required for its observation. 

For a 2D lattice with Vq between 20 and 25Efi (and A'o ~ 170 as in [6]), the conditions 
J ^ U :$> J/Nq are rather well satisfied and the BKT critical temperature, Tbkt ~ J/ks, 
is between 10 and 30nK. Since the BKT transition for the XY model occurs at a tem- 
perature Tbkt ~ J/ks, it should be possible to evidence the transition with measurements 
performed at different temperatures. Using a coarse-graining approach to determine the 
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finite-temperature phase-boundary fine of the Bose- Hubbard model [22] , yields, in 2D and 
for A^o ^ 1 and J/U ^ 1, a critical temperature Tbkt ~ J/ks, in reasonable agreement 
with the XY estimates of Tbkt', in addition, in the chosen range of parameters, the con- 
dition UNq/J — UNq/K ^ 1, satisfies the finite-temperature stability criterion recently 
derived by Tsuchiya and Griffin [23]. 

We notice that the XY Hamiltonian, in the limits A^o ^ 1 and J ^ U ^ J/^qi can be 
also retrieved extending the Gross-Pitaevskii Hamiltonian to finite temperature T » U/ks- 
Replacing the tight- binding ansatz ^{f, t) = Y.ji'jit)^jir) (where ijjj are classical fields and 
not operators) in the Gross-Pitaevskii Hamiltonian one gets the lattice Hamiltonian 

H^-K J:{^I^j + c.c.) + ^^Nj{Nj-1), (9) 
<«J> J 

which is the classical version of the Bose-Hubbard model {Nj =\ ipj p). Writing Nj = 
No + 5Nj, one may neglect the quadratic terms [i.e. (SNj)"^] in the hopping part of the 
Hamiltonian (9), which then reduces in the considered limit to Hxy (8). 

Accurate Monte Carlo simulations yield - for the XY model - the BKT critical tempera- 
ture Tbkt — 0.898 J/Zc^ [24]. When [/ <^ J, a BKT transition still occurs at a slightly lower 
critical temperature Tbkt{U). Intuitively speaking, when U increases, the superfiuid region 
in the phase diagram decreases and, thus, one has Tbkt{U) < Tbkt [20]. An estimate of 
the dependence of Tbkt{U) on U obtained using a renormalization group approach has been 
reported in [25]: when [/ ^ 0, but still U/J <S 1, the only effect of the interacting term 
amounts to renormalize J and to lower the BKT critical temperature [26]. For this reason 
in the following we shall present results from Monte Carlo simulations of the classical XY 
model only. 

In numerical simulations of the finite XY lattice one defines the critical temperature 
as the temperature Tc{L) at which the correlation length equals the size L of the square 
lattice [15]; of course, as L ^ oo, Tc{L) Tbkt- For T = Tc(L), the BKT transition 
is signaled in a finite system by the fact that a suitably defined magnetization, defined by 
Eq.(ll), drops to zero [15]. 
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The emerging physical picture is the following: There are two relevant temperatures for 
the system, the temperature Tbec, at which in each well there is a condensate, and the 
temperature Tbkt at which the condensates phases start to coherently point in the same 
direction. Of course, in order to have well defined condensates phases one should have 

1 /3 

Tbkt < Tbec- The critical temperature Tbec is given by Tbec ~ QMhN^'^uj/kB where 
u) = {CjIuJzY/^ [27]. With the numerical values given in [6], Tbec turns out to be ^ SOOn/T 
for s ;C 20. When T < Tbec, the atoms in the well j of the 2D optical lattice may be 
described by a macroscopic wavefunction ipj. Furthermore, when the fluctuations around 
the average number of atoms per site A^o ^ 1 are strongly suppressed, one may put, apart 
from the factor ^/NQ constant across the array, ipj — e*^?. The temperature Tbkt is of order 
of J/kB' with the experimental parameters of [6] and Vq between 20 and 25Eji, one has that 
Tbkt is between 10 and 30nK, which is sensibly smaller than the condensation temperature 
Tbec of the single well. 

A similar scenario describes also the phases of planar arrays of superconducting Josephson 
junctions [19,20]: they exhibit a temperature Tbcs at which the metallic grains placed on 
the sites of the array become (separately) superconducting and the Cooper pairs may be 
described by macroscopic wavef unctions. At a temperature Tbkt < Tbcs, the array undergo 
a BKT transition and the system - as a whole - becomes superconducting. 

III. OBSERVABLE SIGNATURE OF THE BKT TRANSITION 

In this Section we shall evidence that the experimental signature of the BKT transition 
in bosonic planar lattices is obtained by measuring, as a function of the temperature, the 
central peak of the interference pattern obtained after turning off the confining potentials. 

— * 

Firstly, we show that the peak of the momentum distribution at /c = is the direct 
analog of the magnetization of a finite size 2D XY magnet. In fact, for the XY magnets, 
the spins can be written as Sj — (cos sin 9j) and the magnetization is defined as M = 
(l/N) ■ { I I ) where (• • •) stands for the thermal average; a spin-wave analysis at low 
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temperatures yields M = [2N)-''bT/^^J [15,14]. With discrete BECs at T = 0, all the phases 
9^ are equal at the equilibrium and the lattice Fourier transform of if^j, 

^.— ^E^je-'^ (10) 

exhibits a peak at /c = (/c is in the first Brillouin zone of the 2D square lattice); the 
magnetization is then 

M-dV^ol). (11) 

The intuitive picture of the BKT transition is the following: at T = 0, all the spins point 
in the same direction. As one can see from Fig.2(A), a single free vortex modifies the 
phase distribution far away from the vortex core and, thus, the modulus square of its lattice 
Fourier transform | V^j^ 1^ ^ minimum at A; = 0. At variance, a vortex-antivortex pair 
[see Fig. 2(B)] modifies the phase distribution only near the center of the pair (in this sense 

— * 

is a local defect) and its lattice Fourier transform has a maximum at /c = 0. Analytical 
expressions for the lattice Fourier transform of a single vortex and a vortex-antivortex 
pair may be worked out in detail and are reported in the Appendix B. 

Upon increasing the temperature, vortices are thermally induced. For T < Tbkt only 
bound vortex pairs are present, and on average the spins continue to point in the same 
direction. When the condensates expand, a large peak (i.e., a magnetization) is observed 

— * 

in the central k — Q momentum component, as shown in Fig.2(C). Rising further the tem- 
perature, due to the increasing number of vortex pairs, the central peak density decreases. 
For T ft! Tbkt, the pairs start to unbind and free vortices begin to appear [see Fig.2(D)], 
determining a sharp decrease around Tbkt of the magnetization. At high temperatures, 
only free vortices are present, leading to a randomization of the phases and to a vanishing 
magnetization. In Fig.3 we plot the intensity of the central peak of the momentum dis- 
tribution (normalized to the value at T = 0) in a 2D lattice as a function of the reduced 
temperature /csT/J, evidencing the sharp decrease of the magnetization around the BKT 
critical temperature. 
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Let us now turn our attention to the interference patterns of the expanding condensates: 
after the bosonic system reaches the equihbrium, one may switch off both the harmonic trap 
and the optical lattice. At this time {t = 0) the momentum distribution ip{p,t = 0) is given 
by the Fourier transform of '^{r,t — Qi) — ^jipj^j{r)- one finds 

V^(p,i = 0) = |.(p-)4 (12) 

~ — * 

where is the Fourier transform of the 3D Wannier functions and hk — {px,Py) is the 
momentum projection on the first Brillouin zone of the 2D optical lattice. Using a gaussian 
for the Wannier functions, $(r) oc e~'^^/^°'^e~^^/^'^^e^^^/^'^^ (with az ^ cr = ax = (Jy, since 
ivr ^ (^z), one has $(p) = x{Px)x{Py)x{Pz) ^ e~^^^^~^Py^^'^^^ , where x{Px) ^ e"^^^^/^'*^ and 
similarly for x{Py) ^^'i xiPz)- Since for s ^ 1 a/d — I/tts-^/^, the momentum distribution 
exhibits well pronounced peaks in the centers of each Brillouin zone: these peaks have 

— * 

different heights and the central k — peak is the largest (see Figs. 2 and 3 of [6]). 

— * 

At t = the amplitude of the k — peak of the momentum distribution is simply given 
by the thermal average of ipQ. By measuring the A; = peak (i.e., ( | ipo p)) at different 
temperatures, one obtains the results plotted in Fig. 2. The figure has been obtained using 
a Monte Carlo simulation of the XY magnet for a 40 x 40 array: we find ksTc ~ 1.07 J. 
In Fig. 2 we also plot the low-temperature spin wave prediction [14] (solid line), as well as 
a fit, first used in [15], valid near Tq and derived form the renormalization group equations 
(dotted line). At times different from t — 0, the density profiles are well reproduced by the 
free expansion of the ideal gas: one obtains ^^(p, t) = x{Pz)^~^^^^''^"^^*^^~^'^^^"^'^^^'f{Px,Py,t), 
i.e. a uniformly accelerating motion along z and a free motion in the x — y plane, with 
'fiPx,Py, t) — x{Px)x{Py)'^k^~^'^'^^^^^^*^'^^^ giving the central and lateral peaks of the momen- 
tum distribution as a function of time for different temperatures. 

An intense experimental work is now focusing on the Bose-Einstein condensation in two 
dimensions: at present a crossover to two-dimensional behaviour has been observed for Na 
[28] and Cs atoms [29]. Our analysis relies on two basic assumptions; namely, the validity 
of the tight-binding approximation for the Bose-Hubbard Hamiltonian and the requirement 
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that the condensate in the optical lattice may be regarded as planar [30]. It is easy to see 
that the first assumption is satisfied if, at T = 0, » (where /j, is the chemical potential), 
and, at finite temperature, hujr ~ ksT. The second assumption is much more restrictive 
since it requires freezing the transverse excitations; for this to happen one should require a 
condition on the transverse trapping frequency cu^. Namely, one should have that, at T = 0, 
hcuz ^ 8K and that, at finite temperature, hcUz ^ ksT (since cUz <^ ujr, the latter condition 
also imphes that hwr ~ ksT). In [6] it is Vq — s • ks ■ O-lS/zX and, for s ^ 20, the tight- 
binding conditions are satisfied since lOHz <; 8K/2Trh; furthermore, if cuz = ^tt ■ IkHz, one 
may safely regard our finite temperature analysis to be valid at least up to T ~ 50nK. We 
notice that the experimental signature for the BKT transition for a continuous (i.e., without 
optical lattice) weakly interacting 2D Bose gas [31] is also given by the central peak of the 
atomic density of the expanding condensates. 



IV. CONCLUDING REMARKS 

Our paper analyzes the finite-temperature phase transitions in 2D lattices of BECs. Our 
study evidences the possibility that Bose- Einstein condensates loaded on a 2D optical lattice 
may exhibit - at finite temperature - a new coherent behaviour in which all the phases of 
the condensates located in each well of the lattice point in the same direction. The finite- 
temperature transition, which is due to the thermal atoms in each well, is mediated by 
vortex defects and may be experimentally detectable by looking at the interference patterns 
of the expanding condensates. Our analysis relies on the approximation that the effect of the 
shallow confining harmonic trap in the x — y plane may be neglected: for a tight confinement, 
we expect that interesting mesoscopic effects will come to play, affecting the BKT transition. 

We observe that for a 3D optical lattice (resulting form the potential Vopt{r) — 
Vo[sin^ {kx) -\- siv? {ky) -\- siv? {kz)]) with Nq ~ 100, the system maps on the 3D XY model, 
which is the prototype of the universality class including the A normal-superfiuid transition 
in ^He [32]. Also in this case, the central peak of the interference pattern tends to zero 
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at the critical temperature, Fig.4. We recall, however, that the 3D transition is a simple 
order-disorder transition (whose critical temperature is given in the thermodynamic limit 
by 2.202J/kB [33]), and is not mediated by the creation of vortices- ant ivortices pairs as in 
the 2D case. A natural question arising from the comparison of Figs. 3 and 4 is how to clearly 
characterize the two transitions: to answer, one may observe that the critical exponent /3 
(defined by M oc (Tc - T)^) is /3 a; 0.23 for the 2D XY universahty class [15] and (3 0.35 
for the 3D XY universality class [34] . 

In conclusion, our analysis strengthens - and extends at finite temperature - the striking 
and deep analogy of bosonic systems with superconducting Josephson junction arrays [1]. 
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APPENDIX A: VARIATIONAL ESTIMATES OF THE COEFFICIENTS K AND U 



In this Appendix we derive a variational estimate of the coefficients K and U entering the 
Bose-Hubbard Hamiltonian (5). The tunneling rate K and the coefficient of the nonlinear 
term U are given by 



K^- dr 



(Al) 



and 



U 



(A2) 



where go — ATrh^a/m and $j(r) is the Wannier wavefunction localized at the j-th well and 
normalized to 1. 

For a 2D optical lattice in the x — y plane one has Ve^t — Vopt + ^m, where the optical 
potential is Kpt = Vo[sin^ [kx] -\- sin^ {ky)] and the magnetic potential is Vm — ^i^ri^"^ + 
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y^) + cu'^z^]. Typical values of the parameters relevant for experiments [6] are u!r,ujz ~ 
27r • 10 — 100 Hz, A = 852 nm, for a total number of sites ~ 3000 and an average number of 
particles per site ~ 170. The condition cuz ^ uJr should be imposed in order to have a 2D 
system. To get a variational estimate for K and U we assume that the Wannier wavefunction 
localized at the j-th well is given by 

In Eq.(A3) the parameters ax, Cy, and have to be determined variationally determined, C 
is a normalization constant given by C = (tt^/^ct^ct^)"^''^ and {xj,yj,Zj) are the coordinates 
of the center of the j-th well. Due to the symmetry of the external potential one has to set 
Cx — cfy = CT' furthermore, since ujr ^ ^^r,^z, one has ^ cr. To fix cr and one recalls 
that the Gross-Pitaevskii energy is [27] 



Sm = J df 



(A4) 



In the tight-binding approximation, ^'(r, t) = Y.j 4'j{t)^j{'^)j one finds that the energy 
contribution Ei of the i-th well to the total energy (A4) is E'j / [ | ■0i P ^(V$j)^ + Kxt I 
■01 P $f + ^ I ■0i 1'^ ^i]) ^-iid, since | ipi pRi A^o, Aq being the average value of particle per 
site, one gets 



Ei K, I dr 



(A5) 



Substituting (A3) in (A5), one gets the following approximate expression for the energy in 
the i-th well 

The optimal values of a and cr^ are determined requiring that dEi/da — and dEi/duz — 0. 
It should be stressed that being a the variational width in the directions x and y, it should 
be less than A/2, the distance between minima of the periodic potential. We notice also that 
for realistic values of the parameters, the variational width cr in the x and y directions has 
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a very weak dependence on the mean field term (and therefore to the number of particles 
in each site): in fact, putting a ~ 0.2iJ,ni and ~ bfim, one finds h'^/2ma^ ~ 10~^^erg ^ 
gfo-/Vo/2(27r)^/^(T^(T2 ~ lQ~^^erg, one sees that the last term in Eq.(A6) contributes very little 
to the determination of a. Since Vok'^ ^ ma;^/2, from Eq.(A6) one then finds that the 
dependence of Ei on a is given by Ei{a) ?a iVo(^^/2m(T^ + Vok'^a'^): from dEi/da — one 
obtains a'^ = /2mk'^Vo, which can be cast in the form cr = Ax, with 

and s = Vq/Er. 

We are now in position to have an estimate for K and U : by using the ansatz (A3) in 
Eqs.(Al)-(A2) we get 



K 1/1 



Er 47rWl6 
and 

U 2mgo 



- e-Vi6x^ - se-Vi6x^ (A8) 



(A9) 



En (27r)V2^V,x^' 

The variational procedure used works well also for ID and 3D optical lattices: in a ID 
optical lattice in the x direction one has V^xt = Vopt + Kn; where Vopt = Vq sin^ (kx) and 
Vjn — '^[ojIx^ + a;^(y^ + 2;^)]. Typical values of the parameters relevant for experiments [2] 
are io^ ^ 2t^ ■ 10 Hz, uj^^2t^- 100 Hz, A = 795 nm, for a total number of sites ~ 100 and an 
average number of particles per site ~ 1000. The variational ansatz is still given by Eq.(A3), 
but now (Ty — (Tz = cr_L and (7x = (7 is the variational width in the laser direction. Proceeding 
as before, one obtains that cr = Ax, with x still given by Eq.(A7). The tunneling rate K is 
given by 



(AlO) 



and the nonlinear coefficient U reads U/Eji = 2mXgo/{27iy^'^h'^a'j_x where x is given by 
Eq.(A7). 
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For a 3D optical lattice one has Vext — Vopt + Vm, where Kpt = VJ)[sin (kx) + sin (ky) + 
sin^ (kz)] and Vm — ycj(x^ +1/^+2;^). Typical experimental values from [7] area; ^ 27r-50 Hz, 
A = 852 nm, for a total number of sites ~ 10^ and an average number of particles per site 
~ 1. The variational ansatz is still given by Eq.(A3), but now — Cy — = cr. One 
obtains that a — Xx, with x given as before by Eq.(A7). The tunneling rate K is given by 



(All) 



Er 47rVV16 2 y ^2 2 

and the nonlinear coefficient U reads U/Er — 2mg'o/(27r)^/^^^Ax^. 

We observe that the tunnehng rates K for ID, 2D and 3D optical lattices, given respec- 
tively by Eqs.(AlO), (A8) and (All), can be compactly written as 

Er-Att^xA^^ 2 f '[2 2 ' ' ^""'"^ 

where j — D — 1. 

APPENDIX B: FOURIER TRANSFORM OF A VORTEX ON A LATTICE 

In this Appendix we report the analytical computation of the lattice Fourier transform 
■0^ of single vortex on a lattice described by ipj — exp (iOj) ; namely, one should compute 

= ^ E ^1 e-^'-^ = ^ E (Bl) 
J J 

In Eq.(Bl), j — {jx,jy) denotes the sites of a square lattice (having N sites) and k — {kx, ky) 
the (quasi) momentum, with and ky valued between — tt and vr. The relationship between 
the lattice Fourier transform (Bl) and the momentum distribution ip{p) in real space is 
provided by Eq.(12). 

To simplify matters, in this Appendix wc set the origin O of the coordinates in the 
center of the central plaquette: setting to 1 the lattice length, the lattice sites are labeled 
by j — (ji/2, J2/2) with ji and j2 positive and negative oc?c? integer numbers. 

The phase distribution of a vortex is depicted in Fig. 2 (A) and is such that 6j equals 
the polar angle (f) with respect to the origin O (e.g., at the site (1/2,1/2) 9j equals 7r/4, 
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at the site (3/2, 1/2) is arctan(l/3), and so on; an antivortex would be characterized by 
—(f)). For future convenience, we shall label the site (ji + 1/2, j2 + 1/2) with the index 
£ = (I I + I j2 |)/2: as a consequence, the four sites (1/2, 1/2), (1/2, -1/2), (-1/2, 1/2), 
and (1/2, 1/2) of the central plaquette will be associated to ^ = 1; the eight sites (3/2, 1/2), 
(1/2, 3/2), ■ ■ -, (—3/2, 1/2) will be associated to i — 2. In other words, the sites with index 
i have chemical distance i — 1 from the four sites of the central plaquette. Of course, the 
number of sites with index i is Ai. If, for simplicity, one considers only sites with i going 
from 1 to a maximum value C, then the total number of sites is 



(B2) 



e=i 



To evaluate the lattice Fourier transform (Bl), one may conveniently separate the sum 



in C sums over sites having the same index £; for instance, for the four sites 



(1/2, 1/2), (1/2, -1/2), (-1/2, 1/2), and (1/2, 1/2) of the central plaquette, one gets 



,(1) kx . ky . (1) . kx ky 
sm 01 cos — sm -f — z cos 0i sm — cos -f 
^^2 2 ^^2 2 



(B3) 



where 0i''^^ —n/A and the sum is restricted only to the sites having £ — 1. A straightforward 
generalization of Eq. (B3) shows that for the sum on the 4£ sites having a fixed £ one has 

(2n - 2m + l)kx . (2m - l)ky 



sin 0^-' cos 



sm 



—I cos 



(2n - 2m + l)kx (2m - l)ky 
sm ^ — cos ^ — 



(B4) 



where 



= arctan 



2m 



2n - 2m+ 1 

is an angle taking values between and 7r/2. 

By using Eq.(B4), the lattice Fourier transform (Bl) may be compactly written as 



(B5) 



i'k 



1 



(B6) 
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If in Eq.(B6) one puts ky — 0, one gets 

^^^^ " — ^ — ' ^ ^ 

by introducing the index m — i—m and rearranging conveniently the partial sums in Eq.(B7) 
one obtains 

V'fc = «m sm (B8) 

m=Q ^ 



with 



-4i ^ 2m + 1 



E 



^(^) nikVi ^(2m+l)2+(2n-2?fi-l)2' 
Eq.(B8) clearly shows that for k^ = ^ one has V'o = 0: of course, with a single vortex, the 
magnetization (11) is zero; furthermore, in the thermodynamic limit L — > oo, — 0, in 
agreement the Mermin- Wagner theorem. We observe Eq.(B6) may be easily applied to the 
analysis of the lattice Fourier transform of a vortex-antivortex pair. 
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FIGURES 



10' 

10 20 , 30 40 50 

FIG. 1. J/U as a function of Vq/Er for ID, 2D and 3D optical lattices {z is the number 
of nearest neighbours, which is respectively 2, 4, and 6). The experimental values are taken, 
respectively, from [2], [6], and [7], with an average number of ^"^Rb atoms Nq = 1000 (ID), 170 
{2D) and 1 (3D). 
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FIG. 2. Lattice Fourier transform | tpj^ |^ (for ky = 0) with: (A) a single lattice vortex; (B) a 
lattice vortex-antivortex pair; (C) T = CSJ/fc^; (D) T = l.lJ/kB- Figures (C) and (D) are Monte 
Carlo snapshots after reaching equilibrium. The lattice is 20 x 20. 
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FIG. 3. Intensity of the central peak of the momentum distribution (normahzed to the value 
at T = 0) as a function of the reduced temperature t = ksT/J in a 2D lattice. Empty circles: 
Monte Carlo simulations; solid line: low-temperature spin wave prediction; dotted line: fit near 
Tc ^ 1.07 J /ks as in [15] (in this figure L = 40). 
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FIG. 4. Monte Carlo results for the intensity of the central peak of the momentum distribution 
(normalized to the value at T = 0) as a function of the reduced temperature t = ksT/J in a 3D 
lattice (with 20^ sites). 
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